Estimation of communication-delays through adaptive synchronization of chaos 
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This paper deals with adaptive synchronization of chaos in the presence of time-varying 
communication-delays. We consider two bidirectionally coupled systems that seek to synchronize 
through a signal that each system sends to the other one and is transmitted with an unknown 
time-varying delay. We show that an appropriate adaptive strategy can be devised that is success- 
ful in dynamically identifying the time-varying delay and in synchronizing the two systems. The 
performance of our strategy with respect to the choice of the initial conditions and the presence of 
noise in the communication channels is tested by using numerical simulations. Another advantage 
of our approach is that in addition to estimating the communication-delay, the adaptive strategy 
could be used to simultaneously identify other parameters, such as e.g., the unknown time-varying 
amplitude of the received signal. 
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I. INTRODUCTION 



Over the last twenty years, synchronization of chaos has attracted much attention from the scientific community 
iHa]. An interesting result is the observation that two identical chaotic systems, starting from different initial con- 



ditions can be synchronized on a stable chaotic time-evolution 7 



llj . Even if the systems are slightly non-identical, 



they may nonetheless converge onto an approximately synchronized chaotic time-evolution 



12 



Stochastic syn- 



chronization for two coupled systems or for an arbitrary network in the presence of time delays has been st udied in 
1, 



21| . Applications of synchronization of chaos include but are not limited to secure communication 



identification 



27 



31 



data assimilation 



32 



331, sensors 



34j , information encoding and transmission 



22 



35 



26|, system 



36| , and mul- 



tiplexing 37"]. In this paper, we show that two chaotic systems can be synchronized in the presence of a time- varying 
unknown communication-delay through an appropriate adaptive strategy; moreover, by using this adaptive strategy, 
each system is able to independently formulate a dynamical estimate of the unknown communication-delay. 
It was proposed in 



that the property of a system of being chaotic can be conveniently exploited in control 
applications. In this paper, we consider an identification problem where an unknown time-varying communication- 
delay between two bidirectionally coupled systems is to be estimated. We show that choosing the dynamics of the 



individual systems to be chaotic is convenient in terms of the identification strategy. Moreover, by using our strategy 
we are able to isochronally synchronize the two systems. 

As a reference application, we consider the problem of identifying a time- varying communication-delay between two 
autonomous moving platforms, though we aim at presenting a general methodology rather than a specific application. 
Moreover, our paper provides new insight into the phenomenon of synchronization of chaos, as we are the first ones to 
formulate an adaptive strategy that successfully addresses synchronization of chaos and identification of time- varying 
communication-delays altogether. 

In Sec. II we present the formulation of our problem and as a reference application, we introduce the problem of 
identifying a time- varying communication-delay between two autonomous moving platforms. In Sec. Ill we formulate 
an adaptive strategy to identify unknown communication-delays in the case of non-autonomous systems, which is 
validated by using numerical simulations. In Sec. IV we show that a decentralized agreement protocol can be used 
together with our adaptive strategy to identify the coupling-delays for the case of autonomous systems. The effects 
of the choice of the initial conditions are investigated in Sec. V. In Sec. VI we present a modified adaptive strategy 
that performs better in the presence of noise in the communication channels. Finally, the conclusions are given in 
Sec. VII. 

II. FORMULATION 

A typical scheme for synchronization of two chaotic systems in the presence of communication-delays is the following 



39|-l42| 



i^{t) = F{x,{t),t) + V,[H{xj{t - n)) - H{xdt - n))], (1) 

i = {1, 2} , j = {3 — i) . F : M" x M+ — > K" describes the dynamics of each uncoupled chaotic system, H : M" — > M is a 
scalar output function, Ti is an 7i-dimensional constant coupling vector, the communication-delay is the time that 
it takes for the signal broadcast by node j to be received by node i. The delays in Eq. ^ are assumed constant 
and known. 

Note that Eqs. ([T]) are a system of delay differential equations. Hence the solution Xi{t), i ~ {1,2}, is determined 
by knowledge of the initial conditions for Xi over the time-interval [—Tmax, 0], where Tmax = max(ri, T2). 

In the particular case in which xi ^ X2 ~ x^, the terms in the square brackets of ^ vanish, and the two systems 



3 



evolve on the synchronous solution, 



Xs{t) = F{Xsit),t). 



(2) 



For a given choice of the functions F and H, stability of the synchronized solution depends on the parameters Fi, F2, 
Ti, and T2. In particular, a condition for stability is that the communication-delays be smaller than the characteristic 
timescale of an uncoupled system ([2]) 



Ti < T^, 



(3) 



z = {l,2}. 

If we assume that the communication-delay from i to j is the same as that from j to i, Eq. ([T]) reduces to 



= F{x^{t),t) + r,[H{xj{t - r)) - H{x,{t - r))]. 



(4) 



i = {1,2}, j = (3 — i), Ti = T2 = T. The assumption that the communication-delays are the same in both directions 
holds true for example for those applications that use line of sight communication. 

In this paper, we choose the function _F in ([2]) to generate synchronous (uncoupled) chaotic dynamics. Chaotic 



22 



43 



44j . Adaptive strategies 



signals have been successfully employed in cryptography and in secure communication 
based on synchronization of chaos have also been proposed to dynamically identify the parameters of unknown systems 
27- 



30[. In such applications, it is possible to exploit the specific properties of chaotic signals. In fact, since a chaotic 
time trace never re pea ts in time, it provides an infinite amount of information that can be used by the identification 
strategy (se e e.g.. 



systems 



30 



45|). Moreover, adaptive strategies have been devised to synchronize self-sustained chaotic 



In this paper, we propose an adaptive strategy that exploits the synchronizability property of chaotic systems to 
estimate a time-varying unknown communication-delay. Our goal is twofold: 

i) synchronizing the two systems; 

ii) dynamically and independently estimating at each system the communication-delay. 

To our knowledge, the use of adaptive strategies to simultaneously achieve synchronization and estimate unknown 
communication-delays has not been previously addressed in the literature. A strategy based on synchronization of 



47| and implemented in 48l- 50 1 . 



of unknown real systems 



48- 


5n| 


27 
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Synchronization of chaos has also been used to identify and predict the dynamics 



30|. Adaptive strategies have been proposed to enhance synchronization of coupled 



51 



dynamical systems by acting on the coupling strengths 

parameters have been adaptively estimated in the presence_of time-delays in 
internal delay of an unknown system has been proposed in 



53| and the topolo gy of interconnections 



55 



54|. Unknown 



A technique to estimate the 



45| . The same problem has been studied in 



60| for the 



case of linear systems. In a recent paper 6l\ , synchronization of chaos in the presence of an unknown communication- 
delay is studied between two unidirectionally coupled systems, whose synchronization is guaranteed by a connection 
to an external master system. In this case, the problem of synchronizing the two systems is separate from that of 
identifying the delay. In this paper, we consider two bidirectionally coupled systems, we do not assume connection 
to an external master system, and we formulate an adaptive strategy that can be used to simultaneously synchronize 
the two systems and identify the delay. 

In what follows we present the adaptive synchronization strategy in terms of a reference application. In particular, we 
consider the problem of identifying a time- varying communication-delay between two autonomous moving platforms. 
We assume that two identical chaotic oscillators are installed at each platform and that these chaotic systems seek 
to synchronize via a signal broadcast from one platform to the other (and viceversa). The signal received at each 
platform is transmitted with a time- varying delay, which depends on the relative distance between the two platforms. 
Information on the communication-delay can be used, e.g., to estimate the relative distance. 

Each platform is characterized by a pair of state variables {pi{t), Xi{t)}, i — {1,2}, where Xi{t) € R" is the state 
of an oscillator installed at platform i, and Pi{t) e M™ is the position of the platform i. We assume that the 
two platforms seek to achieve mutual synchronization of the x's dynamics (the dynamics of the oscillators), while 
independently moving along the trajectories Pi{t). In order to achieve synchronization, each platform broadcasts a 
signal which encodes information on the evolution of the x's dynamics. The signals are received with a non-negligible 
communication-delay, which depends on the relative motion of each platform with respect to each other. 

We consider that each uncoupled oscillator is described by Xi{t) — F{xi{t),t). Moreover, we assume that Pi{t) 
typically evolves on a timescale Tp which is much longer than the timescale on which an uncoupled system evolves, 
that is. 



(5) 



5 

This may be achieved, for example, by choosing the individual oscillators to have very fast dynamics, i.e., much faster 
than the dynamics oipi{t), depending also on the specific application of interest. System i = {1,2} broadcasts a signal 
H{xi{t)), which is received by J = (3 — i) with a delay Ti{t), which is a function of both Pi{t) and Pj{t). Under our 
assumption that Tp ^ T^;, we can approximate Ti(t),T2(i) with the same delay, say T(i), with r(t) being a function 
of the distance \\pi{t) —pj{t)\\. 

When the platforms are coupled, the equations for the oscillators become, 

x,{t) =F{x,{t),t) + r,[n{t) - H{x,{t - T^{t)))], (6a) 
r,{t)^H{x,{t-rit))), (6b) 



i = {1, 2}, j = (3 — i). Here, F, H, and Ti are the same as in Eq. (jj]). The received signal at node i, ri{t), propagates 
from node j to node i with a communication-delay T(t). Note that in Eqs. (jS^), we subtract from the received 
signal H{xj{t — T{t))) the internal signal H{xi{t — Tj'(f))), i = {1, 2} , where T^{t) is an estimate of T{t) at node i. 
In the particular case in which xi — X2 ^ Xs and ~ T2 — t, the terms in the square brackets of ([5^) vanish, 
and the two systems evolve on the synchronous solution Note that in this case, integration of the system of 
equations ^ requires knowledge at each time t of the state variables Xi over the time- interval [t — Tmaxit),t\, where 
Tmaxit) = max(r(<),T{(t),T;^(i)). 

A fundamental issue that typically arises when considering communication between two or more autonomous moving 
platforms is that of synchronization of the time-clocks at each platform 

synchronization of non-autonomous chaotic systems has been studied in 64|. In Sec. Ill we will assume the presence 



62 



63[. For example, the problem of 



of two perfectly synchronized internal clocks at the two platforms. This is reflected by the choice of non-autonomous 
chaotic oscillators at the two platforms, Xi{t) — F{xi{t),t), for which the dynamics at the two oscillators explicitly 
depends on the same variable t. As we will see, this requirement is important for the successfulness of the adaptive 
strategy, i.e., for simultaneously synchronizing the two systems and correctly estimating the unknown communication- 
delays. However, in Sec. IV we show that a decentralized agreement protocol can be used to estimate the delay and 
to synchronize the two oscillators for the case that they are autonomous systems. 
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A. Uniqueness of the solution for both the cases of autonomous and non-autonomous oscillators 

We now consider that the uncoupled dynamics at each platform is autonomous, i.e., we replace F[xi{t), t) — ^ F{xi{t)) 
in Eq. (j6^). We show that when such a modification to our scheme is done, we become unable to univocally determine 
the unknown delay r, by requiring that the two systems synchronize, i.e., that H{xj{t — t)) = H{xi{t — t-)) in Eqs. 



A chaotic time trace never repeats in time. Since H{xi{t)) in ([6]) is chaotic, i — {1,2}, the only possibility for the 
received signal H{xj{t)) to cancel out with the internal signal H{xi{t)) over time is that the two chaotic time traces 
are the same one, possibly translated by a fixed time lag equal to A, i.e., 



An example of A-lag synchronization between two chaotic time series (Eq. ^) is shown in Fig. 1. Our goal in 
this paper will be setting A to 0. Fig. 2 shows an example of lag synchronization between two identical P-periodic 
time series. As can be seen, for this case, lag-synchronization is possible for value of the lags equal to A + NP, 
with N being any integer number. Thus using chaotic dynamics at the individual nodes restricts the possibility of 
lag-synchronization to only one lag A. 

By setting the coupling terms in the square brackets of Eqs. ([5]) to zero, and using ([7]), we obtain. 



m 



{1,2}, J -(3-*). 



xi{t) =x2{t- A). 



(7) 



X2{t~T) =X2(t-T{ - A) 



(8a) 



X2{t-T-A) ^X2{t-T!^). 



(8b) 



Again, by observing that X2{t) is chaotic, Eqs. ([5]) reduce to. 



r = t{ + A, 



(9a) 



T + A = r^. 



(9b) 



In Eqs. ([HI, T is a time- varying external parameter. Then we see from ^ that for each value of r, the two systems 
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FIG. 1: The figure shows an example of A-lag synchronization between two chaotic time series generated by integrating Eqs. 
(|19|l . Our goal in this paper we will be setting A to 0. 



may show A-lag synchronization if the following conditions are satisfied, 



A, 



2A = T^-r{. 



(10a) 
(10b) 



A case of interest is that equations (ITOl) are satisfied for, 



(11) 



corresponding to A = 0. Condition dTTI) corresponds to a particular (desirable) solution for (flOl) . for which both the 
estimates r{ and converge on the true value of r, with the two systems synchronizing with lag equal zero. 

In what follows, we will propose an adaptive strategy to evolve the parameters t- in such a way to minimize the 
coupling terms in the square brackets of Eqs. However, we should note that by setting to zero the coupling terms 
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FIG. 2: The figure shows an example of lag synchronization between two P-periodic time series. As can be seen from the 
figure, lag synchronization can be observed for lags equal to A + NP, with A'^ being any integer number. 



in the square brackets of Eqs. a general solution of the type (jlOp will arise, and not necessarily of type (fTTj) . 

Thus to avoid ambiguities in resolving the coupling delays, hereafter we consider two alternative possible remedies, 
one of which applies to nonautonomous systems and the other one applies to autonomous systems. As a first solution, 
we assume that the dynamics of the uncoupled x's in Eqs. ^ are non-autonomous, i.e., they explicitly depend on 
the time variable t. Then, if such dependence on t is appropriately chosen, under the assumption of chaotic dynamics 
for the x's, we can rule out the possibility that Eq. ([7]) is satisfied for A 7^ 0. To better explain our point, let us recur 
to an example. Consider a simple periodic dependence on time of the type, Xi{t) — F{xi{t), sin(27rAi)); then the only 
way Eqs. ([7]) can be satisfied is for A being an integer multiple of A""'^, i.e., A — {0, A^^, 2A^^, ...}. However, not 
all of these solutions will be stable, with stability being typically bounded by a maximum allowed lag A, say A^ax 
(which depends on the autocorrelation time for an uncoupled system, see 4^)- Therefore, by choosing A^^ > A„iax, 



we can ensure that the only possible stable solution of type ([7]) is for A = 0. The advantage of this approach is that it 
uses the received signal ^jp) as the only available information at each platform and from that it is able to reconstruct 
the communication-delay with which the signal propagates from the receiver to the sender. The disadvantage is that 



9 

it requires the presence of two perfectly synchronized internal clocks at the two platforms. 

A second possible solution ensuring A = in Eqs. ([T0|) is to set T2 = r{. In fact, with this condition, Eqs. (fTO|) 
yields Eq. (1111) . However, we note that combining information available at the two different platforms to satisfy 
the condition that be equal to t( may be unpractical. In Sec. IV, we introduce a decentralized agreement protocol 
with this specific purpose. The advantage of this alternative strategy is that it allows to work with autonomous 
dynamical systems at the two platforms, rather than non-autonomous systems. The disadvantage is that it requires 
more information exchanged between the two platforms. 

III. ADAPTIVE STRATEGY 

Our goal is independently estimating at each platform i the unknown communication-delay T{t) from sole knowledge 
of the received signal ri(t), i = 1,2. To this aim, we present an adaptive strategy to evolve Tj'(i), the estimate at 
platform i, to match the unknown communication-delay T(t). At each node i, we consider a potential, 

[H{x,{t - rm - H{xdt - r^mr- (12) 

From the considerations in Sec. IIA, we see that = if (i) xj = Xi and (ii) t[ ~ t, i — 1,2. Therefore, we seek to 
evolve our estimates t^' to minimize 4*^, through the following gradient descent relations, 

r\{t) ^~a,—^= -2a,[H{xj{t - T{t))) - H{xS - Tl{t)))]DH{x,{t - Tl{t)))x,{t ~ rlit)), (13) 

at > 0, DH represents the derivative of the function H with respect to its argument. Note that (|13p involves 
knowledge of the received signal ri(t) = H{xj{t — T{t))), the state Xi at the past time {t — t[), and its derivative Xi 
calculated at time {t — r/). 

We note here that our proposed strategy, Eqs. (j6ll3p . involves the solution of non- autonomous delay differential 
equations with state-dependent delays. In this paper, the stability of our strategy is tested via numerical simulations 
and provides interesting insights on the viability of our proposed approach. Namely, we show its effectiveness for 
different systems, given that the initial guess on the time delay is sufficiently close to the actual initial delay. Further- 
more, the sensitivity to the choice of the initial conditions and to the presence of noise in the communication channels 
is investigated in Sees. V and VI. 
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FIG. 3: We integrate the set of equations (|6I13I14|) with r evolving according to (fT6)) . T = 1, = 0.5, lj = 10 ^. We 
set ai = a2 = 10"''/2, 7i = 72 = 0.15. From top to bottom: time evolution of xii{t) and X2i{t); synchronization error 
— X2i{t)\; time evolution of T(t) (thin black curve) r{(f) (thick dashed grey curve), and T2{t) (thick dotted black curve); 
time evolution of the delay estimation errors, |r(f) — Ti(f)|, |r(t) — r2{t)\. We specify the initial conditions for ((6]) to be constant 
over the time-interval [—fmax,0], where fmax = maxt Tmax{t), and equal to randomly chosen points on the forced Van der Pol 
chaotic attractor. Moreover, we set r{(0) = 1.8, T2(0) = 0.20. 

A. Numerical Experiments 

In what follows, we consider a numerical experiment, for which we choose the dynamics of each uncoupled system 
ii = F{xi{t),t) to be described by the forced Van der Pol oscillator. We restrict the parameters of the coupled 
equations to be such that with no adaption being performed and T[{t) — T2{t) — r(t), the synchronous solution ([2|) is 
stable. This also corresponds to setting a maximum value on the communication-delay T{t), as stability requires that 
T does not exceed the autocorrelation time of an uncoupled system (for more details, see [4^). 

We set Xi{t) = [xa{t),Xi2{t)], 



F{x^i{t),x,2{t),t) = 



Xi2{t) 

-xn{t) + 8.53x,2(i)[l - xn{t)^] + f{t) 



(14) 
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where, 



f{t) — acos(27r/xt), 



(15) 



a = 1.2, /i — 0.1, and we set H{xi{t)) 



We integrate the set of equations (|6ll3p . with 



r(i) = T + CsinM), 



(16) 



uj = 10""^ and we set T = 1, ^ 0.5. We set the initial conditions for ^ to be constant over the time-interval 
[—Tmax,0], where fmax — niaxt Tmaa; (t) and equal to randomly chosen points on the forced Van der Pol chaotic 
attractor. We set t{ (0) = 1.8, r2(0) = 0.20 (note that integration of Eqs. (|6ll3p only requires knowledge of the delays 
T- at the current time t). 

The results of our computations are shown in Fig. 3. As can be seen, the adaptive strategy is successful in 
synchronizing the oscillators' states xi and X2 on a zero-lag synchronous solution, with T[{t),T2(t) converging on the 
true evolution T{t). 

IV. A DECENTRALIZED AGREEMENT PROTOCOL FOR AUTONOMOUS SYSTEMS 

In this section we focus on the case that the two coupled oscillators are autonomous systems, i.e., we replace 
Xi(t) = F{xi{t),t) in Eqs. ([6]) by Xi{t) — F{xi{t)). For this case, the existence of only one solution of type pT|) can 
be ensured by satisfaction of the extra condition that t[ = T2 (see the discussion in Sec. IIA). In order to satisfy this 
additional requirement, we introduce a decentralized agreement protocol between the two platforms that consists in 
each oscillator i communicating to oscillator j = (3 — i) its estimate t^' of the unknown delay. With this modification, 
each platform j receives at each time two signals, that is, the signal rj = H{xi(t — T(t))) and the delay estimate 
T/(i — T{t)). With this extra piece of information available, we reformulate the potential in as follows, 



i.,(t) = ^,{t) + K\T'(t - Tit)) ~ T[(t - T[m\ 



(17) 



i = {1, 2}, j = (3 — z), where k > is an appropriate scalar. Note that in defining P7|) . we have taken into account 
the transmission-time r needed for the estimate rj to propagate from platform j to platform i. We note that, as we 
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see from the relations pi5|) . r evolves on a timescale that is much longer than r itself, thus we do not expect it to 
vary much over the propagation-time of the signal. We observe from (|17p that the potential > by definition. 
Then, in order for the potential (|T7l) to be zero, it is simultaneously needed that equals zero and Tj(t — r) equals 
Tj'(t — T-). Thus we introduce the following gradient descent relation for evolving the estimates r^'. 



(18) 



ai > 0. Note that (|T8)) involves knowledge of the two received signals H{xj{t — T{t))) and rj(t — T{t)), along with 
knowledge of the state Xi, of its derivative Xi, and of the estimate r^' calculated at the past time {t ~ t[). The results 
of our computations are shown in Fig. 4. We choose the dynamics of each individual system to be described by the 
Rossler autonomous equation, Xi{t) = [xii{t),Xi2{t),Xi3{t)], 



F{xii{t),Xi2{t),Xi3{t),t) = 



-xa{t) - xa{t) 
x,i{t)+{).2xa{t) 

X^5{t)[x^l{t) - 7] + 0.2 



(19) 



We integrate the set of equations (|6ll9ll8p with H{xi{t)) = Xii{t), Ti — [7^,6,0]. The unknown delay r evolves 
according to Eq. ^ with T = 0.3, C = 0.15, w = lO^^. We set ai = aa = 10-'*/2, 71 = 0.6, 72 = 0.5, and k = 1. We 
select the initial conditions for xn, Xi2, xa to be constant over the time-interval [—fmax, 0], where fmax = niax( Tmaxit), 
and equal to randomly chosen points on the Rossler chaotic attractor, t((0) ~ 0.65, T2(0) ~ 0.10. As can be seen, 
the adaptive strategy along with the decentralized agreement protocol is successful in synchronizing the oscillators' 
states xi and X2 on a zero-lag synchronous solution, with T[{t) and T2{t) converging on the true evolution T{t). 

Furthermore, we have tested the performance of the adaptive strategy described in this section with respect to 
different sets of parameters. We have run numerical simulations in which we integrate the set of equations (|6I18I19P 
with K — 1 and r evolving according to (IT5|) . T = 0.4, ^ = 0.2, to = 10^^. In order to assess the effectiveness of the 
adaptive strategy we monitor the estimation error. 
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FIG. 4: We integrate the set of equations (|6I18I19|I with r evolving according to ([TeJ, T = 0.3, C = 0.15, a; = 10"^. We set 
a\ = 0L2 — 10~*/2, 7i — 0.6, 72 = 0.5, k = \. From top to bottom: time evolution of x\\(i^ and a:2i(t); synchronization error 
\x\\{t') — X2i(t)j; time evolution of rii) (thin black curve) T{(f) (thick dashed grey curve), and T2(t) (thick dotted black curve); 
time evolution of the delay estimation errors, |r(f) — Ti(t)|, |r(f) — ''"2(i)|. We select the initial conditions for Xi\,Xi-i,Xiz to be 
constant over the time-interval [— r„iai,0], where fj^ax = maxt Tmaa;(t), and equal to randomly chosen points on the Rossler 
chaotic attractor, r{(0) = 0.65, r2(0) = 0.10. 

In Fig. 5 we plot the average estimation error. 



/.2X10'' 

1-4 



>=10-M E,{1i)di, (21) 

I = 1, as we vary the coupling gains 71 — 72 between the two systems and the adaptation gains a\ = a2- As can be 
seen, the strategy is quite sensitive to the choice of both the coupling gains and the adaptation gains, which suggests 
that these parameters should be carefully tuned in order to obtain a desired or improved performance. 

We observe that the use of the approach proposed in this section is not limited to the reference application of 
identifying a time-varying communication-delay between two autonomous moving platforms. Indeed, the version of 
the adaptive strategy that applies to autonomous systems does not require the presence of synchronized clocks at the 
two platforms; hence, it could be used to synchronize two autonomous clocks in the presence of clock drifts or skews 
(see 65| for a survey on clock synchronization) or to estimate propagation delays in an asynchronous code-division 
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FIG. 5: We integrate the set of equations (|6I18I19|I with r evolving according to (fT6)) with T = 0.4, ( = 0.2, u = 10~^. We 
set k = 1. (a) We plot the average estimation error < _Ei > as we vary the coupling gains between the two systems 71 — 72 
with ai = a2 = 10~''/2. (b) We plot the average estimation error < _Ei > as we vary the adaptation gains ai = 02 with 
71 = 72 = 0.6. In both plots, each point is an average over many different realizations (i.e., over different choices of the initial 
conditions for xi and 2:2). 



multiple access communication system 66|. As another possible application, our setup could be used as a bistatic 
radar, where one platforms (acting as a transmitter) communicate a signal to the other one (acting as a receiver) in 
the presence of a moving target. In this case, knowledge of the time delay for the signal to bounce off the moving 



target and reach the receiver could be used to extract information on the position of the target. 



V. EFFECTS OF THE INITIAL CONDITIONS 

In this section, we investigate the effect of the choice of the initial conditions {r((0), T2(0)} on the successfulness of 
our adaptive strategy. We consider both the forced Van der Pol equation ([14)) as well as the forced Rossler equation 
(IT^ to describe the dynamics of an uncoupled system. We choose t((0) = t(0) + Ai, T2(0) — t(0) + A2, where is 
a measure of how off our initial guesses r^' are with respect to the true value of t(0). We select the initial conditions 
for Xi to be constant over the time-interval [—fmax,0], where fmax = maxt Tmaa: (i) , and equal to randomly chosen 
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points on the Rossler chaotic attractor. We run numerical simulations in which equations (|6I13I14I) are integrated for 
a long time, H{xi(t)) — Xi2{t), Ti — [0,7i], ai = a2 — 10~*/2, and 71 = 72 = 0.15. At each time t, we monitor 
the estimation error < E2 > defined in Eq. ((2T|) versus the initial mismatch A2 for different values of Ai: circles 
correspond to Ai = —0.9, squares correspond to Ai = 0, and asterisks correspond to Ai — 0.9. Each plotted point is 
averaged over many different realizations, i.e., over different choices of the initial conditions for xi and X2 (constant 
over the time-interval {[—fmax, 0]) ). As can be seen, for < A2 ^ 1.2, our strategy is able to produce good estimates 
of the unknown time- varying delay T{t) within the time-span of our simulations (2 x 10^). We remark that for Ai 
in the interval —0.5 < Ai < 0.5 (which represents a variation of up to ±50% on our initial guess with respect to the 
true value of r(0) — 1), the results of our computations are very similar to the case of Ai = 0. 

In the inset we investigate the effects of the initial conditions on the decentralized agreement protocol introduced in 
Sec. IV. We take the individual systems to be Rossler oscillators (|19p : we set Ai = 0.8 and plot the average estimation 
error < E2 > versus the initial mismatch A2. As can be seen, for —0.4 < A2 ^ 1.5, the adaptive strategy is able to 
produce good estimates of the unknown time- varying delay T{t) within the time-span of our simulations (2 x 10*). 



In this section we consider the presence of noise in the communication channels between the two platforms. We 
assume additive noise and we replace in ^p) by. 



where ei{t) is a zero- mean independent random number of unit variance drawn from a Gaussian distribution, ai is a 
multiplicative factor, and s is the time step of our integration method. 

To better deal with the presence of noise in the received signals, we propose an alternative formulation of our 
adaptive strategy for which the potential p2)) is replaced by. 



From (|23p we see that ^'^ is an exponential moving average of the squared synchronization error with averaging time 
i/^^. We require to be larger than T^, the characteristic timescale on which an uncoupled system evolves, and to 



VI. EFFECTS OF NOISE 



r,(t) = H{xjit - Tit))) + V^aMt), 



(22) 




(23) 
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FIG. 6: We integrate the set of equations (|6I13I14|I with r evolving according to Eq. with T = 1, C = 0.5, u = 10"'^. The 
main plot shows the average estimation error < E2 > versus the initial mismatch A2. Diflterent markers are for different values 
of Ai: circles correspond to Ai = —0.9, squares correspond to Ai = 0, and asterisks correspond to Ai — 0.9. Each point is an 
average over many different realizations (i.e., over different choices of the initial conditions for xi and X2); ai = 02 = 10~^/2, 
71 = 72 = 0.15. The inset shows the results of simulations that study the effects of initial conditions on the decentralized 
agreement protocol introduced in Sec. IV. We integrate the set of equations (|6I18I19P with r evolving according to (|16[) . T — 0.4, 
— 0.2, uj — 10~^ . We set ai = 02 = 10~''/2, 71 — 0.6, 72 = 0.5, k = 1. The inset shows the average estimation error < E2 > 
versus the initial mismatch A2 for Ai = 0.8. 



be smaller than Tp, the timescale on which the communication-delay changes, 



< i^-i < Tp. (24) 



Note that and = only if xi = X2 and = r (see Sec. IIA for a more detailed discussion). 
Follovifing our previous derivations, the equations for the alternative adaptive strategy are, 



r'z = -a,— 4 =-2a,^„ (25) 
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where > 0, and 

^^{t)^l^ J e-'^'^'-'^niO) - H{x,{e - T[{e)))]DHx,{e - Tm)dO. (26) 

We observe that the quantity in (|26)) obeys the following differential equation, 

Ut) = -<^it) + '^hW - H{x^it - Tl{t)))]DHx^{t - Tl{t)). (27) 

Thus our alternative adaptive strategy is completely described by the set of Eqs. ((6k l22l25l27l) . We note that with 
respect to the original formulation, the alternative strategy requires integration of an additional differential equation 
at each platform. 



0.8 




FIG. 7: The figure shows the results of numerical integration from t — Q to t = 10* of the set of equations (|6ll9l22l28p . 
corresponding to our alternative formulation of the adaptive strategy for autonomous systems, H{xi{t)) = Xii{t) and Ti — 
[7i,0,0]. r{t) evolves according to Eq. with T = 0.3, ( = 0.15, uj = 10"^ We set ai = Q2 = 10"*/2, 71 = 0.6, 72 = 0.5, 

K = 1. We record the average estimation error < E2 > for increasing values of the noise ratio o /a^, where a = g\ = a2 and 
<Jx = 5.5 is the standard deviation of the time evolution of H{xsi). Squares correspond to v = 1.5 and diamonds to 1/ — 3. 
Circles correspond to runs in which the same simulation is repeated for our original formulation of the adaptive strategy, Eqs. 
(|6I13|) . Each point is an average over many different realizations, i.e., over different choices of the initial conditions for xi and 
X2 (constant over the time-interval [— rmai,0]). 
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We have also considered the effects of noise on the decentrahzed agreement protocol introduced in Sec. IV. For this 
case, Eq. (fT8|) is replaced by 



(28) 



Vt^ + {[n{t) - H{x,{t - T[{t)))]DH{x,{t - Tl{t)))x,{t - Tl{t)) 



j = 1, 2, J = 3 — i, Q!i > 0, Tj{t — T{t)) — t' j(t — T{t)) + ^/s(Jiei{t), with s, cr,;, and being the same as in Eq. (j22p . 

In Fig. 7 we show the results of numerical experiments involving two coupled chaotic Rossler oscillators. We 
integrate the set of equations (|6ll9l22l28p . corresponding to our alternative formulation of the adaptive strategy for 
autonomous systems. We consider that the communication-delay r slowly varies in time according to Eq. (|16l) . with 
T = 0.3, C = 0.15, uj = 10^"^ (see the figure caption for the simulation parameters). We record the average estimation 
error < E2 > for increasing values of the noise ratio a/a^, where a ^ ai ^ a2 and — 5.5 is the standard deviation 
of the time evolution of H{xsi)- Squares correspond io v — 1.5 and diamonds to = 3. Circles correspond to runs 
in which the same simulation is repeated for our original formulation of the adaptive strategy, Eqs. (|6ll8p . As can be 
seen, for low noise (low values of cr) the original adaptive strategy and the modified adaptive strategy (for both cases 
of = 1.5 and v = i) produce similar results. As a is increased, the average estimation error < E2 > grows. For 
large enough noise <t {a }^ 0.12(T^), the modified adaptive strategy outperforms the original adaptive strategy. 

VII. A STRATEGY FOR SIMULTANEOUSLY ESTIMATING THE COMMUNICATION-DELAY AND 

THE SIGNAL AMPLITUDE 

In the previous sections we have shown that a suitable adaptive strategy based on synchronization of chaos can 
be used to estimate an unknown communication-delay between two coupled systems and achieve synchronization 
between them. As a reference application, we have proposed that such a strategy could be used to estimate the 
relative distance between two autonomously moving platforms. 

Another aspect that we have not taken into account is that since communication occurs between two moving 
platforms, the amplitude of the signal received at each platform could be unknown and time-varying as well. One 
viable strategy could be to try to infer the unknown amplitude from the delay. In this section, we will show that the 
adaptive strategy proposed in Sees. IV can be effectively extended to simultaneously and independently estimate the 
unknown communication-delay and the unknown amplitude, even when they evolve independent of each other and in 
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a noisy environment. We wish to emphasize that an adaptive str ategy to estimate the amphtude of the signal alone 



has already been studied in 
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48|-l50|. 



47| and experimentally tested in 
In order to take into account the effect of the unknown time-varying amplitude of the signal received at each 
platform, we replace equation ^ by the following, 



x,{t) ^F{x,{t)) + T,[n{t) - A[{t)H{xdt - Tl{t)))l (29a) 
n{t)^A{t)H{xj{t-T{t))), (29b) 

i = {1,2}, j = (3 — i), where both the amplitude A{t) and the communication-delay T{t) are unknown and slowly- 
evolving in time and A'^{t) is an estimate at platform i of the unknown amplitude A{t). Again, we assume that the 
two platforms are coupled in line of sight, so that it can be assumed that the attenuation of the signal from platform 
i to platform j is the same as that from platform j to platform i and that the communication-delay is the same in 
both directions. 

We now seek to obtain an adaptive strategy to simultaneously and independently estimate at each platform both 
the unknown amplitude A(t) and the communication delay T{t). We assume that the two chaotic systems at each 
platform arc autonomous, thus we use a decentralized agreement protocol similar to the one presented in Sec. IV. 
Similarly to Sec. IV, we introduce at each platform i a potential, 

*.(i) = [A{t)H{x,{t - T{t))) A[{t)H{xdt - TimY + K[A{ty^{t T{t)) - A[{t)T[{t - Ti{t))Y- (30) 

Note that > Q and '^i = only if xi = X2, t[ = r, and A'^ = A. The latter follows from the observation that 
H(xi) and H{x.j) are chaotic and rapidly evolving in time, while A{t) and A[{t) can be considered constant over the 
timescale on which chaos evolves. Then the only way that the first term on the right hand-side of pOI) can be set 
equal to zero independent of time is for A'^ — A, for which case A[ = A can be factored out of the square brackets in 
(1301) . Then the condition for ^E'i = is that both the terms in the square brackets are equal zero, yielding xi = X2, 
and Tj' = T (for more details, see Sec. IIA). 

Thus we seek to formulate an adaptive strategy that seeks to evolve A[{t) and Tj'(f) in such a way to minimize the 



20 



potential (pO)) . To this aim, we introduce the following gradient descent relations, 



A[{t)H{x.{t-Tm)]H{x.{t-Tm) 



'i 



+ n[A{t)T'^{t-T{t)) 



A',(t)Ti(t~Timriit-Tim, 



(31a) 




2a,{[A{t)H{x,{t-T{t))) 



A[{t)H{x,{t ~ T[m]A[{t)DH{x,{t - rim^^it " r'm 



K\A{t)Tr{t~T(t)) 



m)r[(t~T[mA!m. 



(31b) 



with tti, /3i > 0, i = {1, 2}. 

We have numerically tested the adaptive strategy described by equations (I29l3ip with T(t) evolving according to 
Eq. with T = 0.3, ( ~ 0.15 and A{t) evolving according to the following equation. 



with A = 1, T] = 0.1, LUA = 5 X 10 ^. We choose each individual system to be a chaotic Rossler oscillator (|19l) . 
H{x,{t)) = xa{t) and = [7»,0,0]. We set ai = = /3i = /32 = lO""*, 71 = 0.6, 72 = 0.5, k = 1. The results 
of our numerical experiment are shown in Fig. 8. As can be seen, after a transient the adaptive strategy is able to 
track the time-evolutions of both T{t) and A{t). Also, the time evolutions of xi{t) and X2{t) are synchronized (not 
shown). Note that the strategies to estimate the unknown delay and amplitude (given in Eqs. (PTb ) and (15Tb )) are 
run simultaneously and independently of each other. 



Our strategies described in this paper take advantage of the properties of chaotic sequences of never repeating 
in time. Another way of generating long non-repeating time-traces is to consider sequences of independent random 
numbers. For example, consider the A^-random sequence 



A{t) = A + r/sin(wAi), 



(32) 



VIII. DISCUSSION 




0, with probability p, 



1, with probability 1 — p. 



(33) 
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FIG. 8: The figure shows the results of numerical integration from t = 0tot = 2x 10* of the set of equations (|29I19I31|I . 
corresponding to the version of the adaptive strategy for simultaneously estimating the communication-delay and the signal 
amplitude, H{xi{t)) = Xii{t) and Ti = [7^,0,0]. r(t) evolves according to Eq. (|16p . with T = 0.3, ^ = 0.15; A{t) evolves 
according to Eq. with A = 1, 77 = 0.1, uja = uj/2 = 5 x 10"''. We set ai = aa = ft = /32 = lO"'*, 71 = 0.6, 72 = 0.5, k = 1. 
From top to bottom: time evolution of r(f) (thin black curve), Ti{t) (thin grey curve), and T2{t) (thick dotted black curve); time 
evolution of the delay estimation errors, |r(f) — r[{t)\, |r(t) — T2{t)\; time evolution of A{t) (thin black curve), A'i{t) (thin grey 
curve), and A'2{t) (thick dotted black curve); time evolution of the amplitude estimation errors, \A{t) — A'i{t)\, \A{t) — A'2{t)\. 
We select the initial conditions for xn, Xi2, Xis to be constant over the time- interval [— Tm,aa:,0], where fmax ~ maxt Tmax{t), 
and equal to randomly chosen points on the Rossler chaotic attractor, t[{0) = 0.7, T2(0) = 0, A'i{0) = 0.6, ^2(0) = 1.4. 

n = 0,1, N , with < p < 1- Then if we select two s-independent random sequences as defined in (j33p . we can 
compute the probability P that they arc the same, P = (1 — 2p{l — p))^. It is easy to see that unless p = or p = 1, 
P tends to as iV tends to infinity. Hence it is very unlikely to produce two random sequences that are the same. 

A way of getting around this problem, is to rely on a pseudorandom number generator, that is, a deterministic 
algorithm for generating a sequence of numbers that approximates the properties of random numbers. Such an 
algorithm requires initialization by a relatively small set of initial values, called the seed. Two sequences that are 
initialized by the same seed are identical. The seed plays the same role as the initial condition of a chaotic dynamical 
system and indeed a way of producing pseudorandom sequences is to use an underlying deterministic chaotic process 
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69( 1 . In this respect, the use of a pseudorandom algorithm could represent a valid alternative to the use of a 



chaotic dynamical system for the delay-identification problem discussed in this paper. 
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However, the key-property that we present and exploit in this paper is that two chaotic dynamical systems, even 
if initialized from two different initial conditions, can be synchronized, i.e., by introducing an appropriate feedback 
mechanism, they can be maintained in a stable synchronized chaotic time-evolution l|,|7|,|8|. Furthermore, it has been 
shown that even in the presence of slight non-identicality between the coupled systems and slight deviations from 
nominal conditions, a stable approximately synchronous evolution can be reached and maintained in time 
Thus the use of a closed- loop architecture has the following advantages: 

(i) The two systems do not need to be initialized with the same initial condition, i.e., a-priori knowledge of a seed 
is not required. 

(ii) The mechanism is robust with respect to noise and slight deviations from nominal conditions, such as slight 
non-identicality in the individual systems, slight variations in the amplitude of the received signal, or slightly different 
environmental conditions. 

To conclude, the advantage of our chaos-synchronization strategy is that it relies on a closed-loop architecture that 
allows to maintain synchronization between two distant systems even in the presence of noise and non-identicality in 
the individual systems or mismatches in the initial conditions. 

On the other hand, noise or small mismatches in the parameters of the individual chaotic systems being coupled 
can be responsible for the onset of bubbling 
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70|-|72||. i.e., rare intermittent large deviations (bursts) from synchro- 
nization. Bubbling occurs when the synchronized state is stable for typical chaotic orbits but is unstable for certain 
unstable periodic orbits within the synchronized chaotic attractor. Consider for example the set of Eqs. (|6ll8p . 
For the case of constant t, these equation allow an invariant set, the so-called synchronous manifold (SM), defined 



as xi{t) = X2{t), t[ 



T. Stability of the SM with respect to infinitesimal transversal perturbations can be 



quantified in terms of the maximum transverse Lyapunov exponent (MTLE) of the system. Now assume the MTLE 
is negative, implying that stability is observed for any initial condition on the attractor. However, in the presence of 
either noise, or small mismatches in the parameters of the individual systems, or even a slowly time- varying r(t), it is 
possible that the trajectory eventually gets close to an unstable periodic orbit embedded in the attractor having an 
associated positive MTLE. If this happens, the trajectory may be repelled away from the SM and eventually return 



7C 



7l|. In 



lai bub- 



close to the SM after some time, giving rise to a burst. This phenomenon is called bubbling 
bling is observed for a network of coupled dynamical systems, each of which independently implements an adaptive 
strategy to maintain synchronization. Analogously, we expect bubbling to eventually arise (depending on the choice 



of the equation parameters), for the problem described by Eqs. (j6ll8p . i.e., for a case in which two coupled systems 
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implement an adaptive strategy to estimate a communication-delay. In particular, bubbling is likely to arise in the 
case in which the unknown delay is time- varying, as mismatches in the delay estimate \T'i{t) — T(i)|, no matter how 
small they are, may generate it. 



IX. CONCLUSIONS 



In this paper, we considered synchronization of two bidirectionally coupled chaotic systems that communicate 
with an unknown time-varying delay. We presented a nonlinear adaptive strategy to simultaneously identify the 
time-varying communication delay and achieve synchronization. 

To achieve the twofold goal of synchronizing the chaotic systems and identifying the unknown delay, we proposed 
an adaptive strategy based on the minimization of an appropriately defined potential through a gradient descent 
technique (see also 
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48|). The proposed strategy is presented in two alternative formulations to deal with 
either non- autonomous or autonomous systems. The version for non-autonomous systems requires the presence of 
two perfectly synchronized internal clocks at the two platforms, while the version for autonomous systems is based 
on a decentralized agreement protocol, i.e., on exchanging information on the respective delay-estimates between the 
two platforms. 

As a reference application, we have considered the problem of identifying the communication-delays between two 
autonomous moving platforms. We have assumed that two identical chaotic oscillators are installed at each platform 
and that these chaotic systems seek to synchronize via a signal broadcast from one platform to the other (and 
viceversa). Moreover, we have assumed that the signal is transmitted with an unknown communication-delay that 
depends on the distance between the platforms. Our strategy could be used to dynamically estimate at each platform 
the distance at which the other platform is at any given time. 

The effectiveness of the proposed approach has been illustrated by means of extensive numerical simulations, 
showing that the adaptive strategy can be effective in synchronizing the chaotic systems and in correctly estimating 
the time-varying communication-delay. Another advantage of our approach is that in addition to estimating the 
communication-delay, the adaptive strategy could be used to simultaneously identify other parameters (see e.g.. 
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47[). For example, in Sec. VII we have shown that the adaptive strategy is effective in simultaneously and 
independently estimating the unknown time- varying communication delay and the unknown time- varying amplitude 
of the received signal. 

To investigate the robustness of the proposed approach, we have numerically tested the performance of our strategy 



24 

with respect to the choice of the initial conditions. We have also investigated the effects of noise in the communica- 
tion channels between the platforms and proposed an appropriately modified adaptive strategy for which a relevant 
enhancement of the performance has been observed in the presence of noise. 

Our adaptive strategy (especially in its version for nonautonomous systems) could be useful in applications as dif- 
ferent as the identification of communication delays between moving platforms, clock synchronization, the estimation 
of propagation delays in an asynchronous code-division multiple access communication system, and the localization 
of a mobile target through a bistatic radar system. However, a main limitation to our approach seems to be that for 
some applications (e.g., clock synchronization in a computer network) the communication delays may be asymmetrical 
(i.e., different in the two directions from i to j and from j to i). The condition of symmetrical delays, which holds 
for line of sight communication, seems to be relevant to the formulation of our strategy that applies to autonomous 
systems. How to extend the strategy for autonomous systems to the case of asymmetrical delays is the subject of 
ongoing investigations. One viable approach would be to increase the amount/type of information that the platforms 
are allowed to exchange. As a reference for future work, we also note the importance of testing our strategy with 
systems displaying larger complexity than the simple Van der Pol and Rossler oscillators considered in this paper 
(e.g., systems of interest for practical applications). 
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